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The oblique impacts of nanoclusters is studied by means of Molecular Dynamics and theoreti- 
cally. In simulations we explore two models - Lennard- Jones clusters and particles with covalently 
bonded atoms. In contrast to the case of macroscopic bodies, the standard definition of the normal 
restitution coefficient yields for this coefficient negative values for oblique collisions of nanoclusters. 
We explain this effect and propose a proper definition of the restitution coefficient which is always 
positive. We develop a theory of an oblique impact based on continuum model of particles. A 
surprisingly good agreement between the macroscopic theory and simulations leads to the conclu- 
sion, that macroscopic concepts of elasticity, bulk viscosity and surface tension remain valid for 
nanoparticles of a few hundreds atoms. 

PACS numbers: 45.70.-n,36.40.-c,45.50.Tn 
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Introduction. Inelastic collisions, where a part of me- 
chanical energy of colliding bodies transforms into heat, 
are common in nature and industry. Avalanches, rapid 
granular flows of sand, powders or cereals may be men- 
tioned as pertinent examples [l|, |2[. Moreover, inelas- 
tic collisions define basic properties of astrophysical ob- 
jects, like planetary rings, dust clouds, etc. An important 
characteristic of such collisions is the so-called normal 
restitution coefficient e. According to a standard defini- 
tion, it is equal to the ratio of the normal component of 
the rebound speed, g' (prime states for the post-collision 
value), and the impact speed, g 
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gn 
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The unit inter-center vector n = ri2/|ri2| at the colli- 
sion instant (ri2 = ri — r2) specifies the impact geome- 
try. Since particles bounce in the direction, opposite to 
that of the impact, e is positive, e > 0, and since the 
energy is lost in collisions, e is smaller than one, that is, 
< e < 1. This is a common statement in the majority 
of mechanical textbook, where it is also claimed that e is 
a material constant. Recent experimental and theoreti- 
cal studies show, however, that the concept of a restitu- 
tion coefficient is more complicated: First, it depends on 
an impact speed [3|-|5[, second, it can exceed unity for a 
special case of oblique collisions with elastoplastic plate 
[6|, where the energy of normal motion can increase at 
the expense of the energy of tangential motion [y]. Still, 
it is believed that e < 1 for a true head-on collision. 
The concept of a restitution coefficient, as a basic one 
of the classical mechanics, has been introduced long ago 
by Newton; it addresses an impact of macroscopic bod- 
ies. The increasing interest to nanoparticles, inspired by 
its industrial significance, raises an important question, 
to what extent the macroscopic concepts are applicable 
and whether they acquire new features at a nanoscale. 
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FIG. 1: (Color online) Initial (left) and final stage (right) 
of the nanocluster collision. The initial relative velocity is 
Vi2(0) = V and the incident angle is 7. The unit normal n 
specifies the orientation of the contact plane. For large 7 a 
noticeable reorientation of this plane is observed. Here the 
collision of H-passivated Si nanospheres (model B) is shown. 



The collisions of nanoclusters has been studied in de- 
tail numerically |7H1C| . It was observed that the surface 
effects, due to the direct inter-cluster van dcr Waals in- 
teractions, play a crucial role: The majority of collisions 
of homogeneous clusters, built of the same atoms, lead to 
a fusion of particles [7| ; they do not fuse for high impact 
speeds, but disintegrate into pieces [7|. This complicates 
the analysis of restitutive collisions, which may be more 
easily performed for particles with a reduced adhesion. 
Among possible examples of such particles are clusters 
of covalently bonded atoms, especially when their surface 
is coated by atom of different sort, as for H- passivated 
Si nanospheres [8|. These particles can rebound from a 
substrate, keeping their form after an impact unaltered 
|8|. The bouncing nanoclusters demonstrate a surpris- 
ing effect - the normal restitution coefficient can exceed 
unity even for strictly head-on collisions [9|. 

In this Letter we investigate the oblique impact of nan- 
oclusters with the reduces adhesion by means of Molecu- 



lar Dynamics (MD) and theoretically, using concepts of 
continuum mechanics. Unexpectedly, we have found that 
the normal restitution coefficient, as defined by Eq. ([T|), 
acquires for large incident angles negative values, e < 0. 
We explain this effect by the reorientation of the contact 
plane during an impact and quantify it. Moreover, we 
propose a modified definition of e, which preserves its 
initial physical meaning and yields always positive val- 
ues. To describe the collision of nanoclusters we develop 
a continuum theory. Surprisingly, the macroscopic ap- 
proach quantitatively agrees with MD even at nanoscale. 
MD simulations. We study two models - a sim- 
plified model (A), which mimics interactions of nan- 
oclusters with the reduced adhesion and realistic model 
(B) for interaction of nanoclusters with covalently 
bonded atoms - H-passivated Si nanospheres. For 
the model A, which is less computationally expensive, 
we adopt the Lennard- Jones ( LJ ) potential (j){r) = 
4e [{cLj/ry^ — (cTLj/r)^] for the interaction between two 
atoms in the same cluster and the modified LJ potential 
<l)\nt{r) = 4e [(o'Lj/'')^^ - c(crLj/?')'^] for the interaction 
between atoms in two different clusters. Here the cohe- 
sive parameter c = 0.2 controls the adhesive force [9|, lid ] 
between clusters, while e, ctlj, and r are, respectively, 
the depth of potential well, the diameter of the repulsive 
core, and the distance between two atoms. The nan- 
oclusters of A^ = 500 atoms were prepared by the two 
step temperature quench to obtain a rigid amorphous 



particle 12| . The diameter of nanocluster d was defined 



as the maximum distance between the center of mass of 
the nanocluster and the atom on the surface, so we find 
d = 10.46 (Tlj- For the model B we adopt the Tersoff 
potential [13[ for the covalent Si-Si, Si-H, and H-H bond- 
ings. The Si nanospheres, containing 2905 Si atoms ar- 
ranged in a diamond structure, are fully coated by 852 H 
atoms. The radius of Si nanosphere is about d = 2.4 nm. 

We fix the modulus of the relative inter-cluster velocity 
vi(0) - V2(0) = vi2(0) = V and set it to T^ = 1.0 ^/e/^ 
and 1850 m/s for the model A and B, respectively. We 
vary the incident angle 7 between n and V (see Fig. 
[T|), so that the normal impact velocity, Vn = ycos7 is 
changed. The nanoclusters do not rotate before an im- 
pact and have zero angular velocities, a;i(0) = a;2(0) = 0. 
To make an ensemble average, we randomly turn one of 
the clusters around the axis, passing through its center 
and perpendicular to the contact plane. Due to rough 
atomic surfaces of the clusters, this results in varying 
contact configurations at each impact. Hence, for every 
incident angle 7 we perform averaging over 100 collisions 
with different contact conditions for model A and over 
10 collisions for model B. The clusters' deformation dur- 
ing an impact is quantified by the normal displacement, 
^„(t) = d— |ri2(t)| = d — ri2{t). Wc define the beginning 
of a collision at t ~ Q and the end at t = t^ through the 
conditions, ^„(0) = ^„(tc) = 0. 

Simulation results for the normal restitution coefficient 
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FIG. 2: Dependence on the incident angle 7 of the normal 
restitution coefficients e, and e according to the standard def- 
initions @ and modified definition ((S)). Open circles and 
squares are respectively the MD results for e and e, while 
dashed and solid lines correspond to theoretical predictions. 
Upper panel refers to the model A and lower panel - to the 
model B. Note that the coefficient e is always positive. 



for the models A and B are shown in Fig. [2] (upper and 
lower panel respectively). As it is seen from the figure, 
the restitution coefficient e, defined by Eq. ([T]) becomes 
negative for large incident angles 7. Such unusual be- 
havior of e at nanoscales may be understood if we notice 
that the orientation of the contact plane, characterized 
by the unit vector n(t) — ri2{t)/ri2{t), significantly al- 
ters during the collision. Fig. [1] This is quantified by the 
angle a between the initial and final orientations of n(t). 



cos a = n(0) • n(tc) 



(2) 



The dependence of a on the incident angle 7 is shown in 
Fig. [31 If a is large, the normal restitution coefficients 
becomes negative. Fig. [21 

Modified definition of e: To analyze this effect, con- 
sider the relative velocity of particles at their contact. 
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g = vi2-h- [n X 1^12] = -^„n-fri2n+- [n x ^12] , (3) 



+ CJ2 and we use V12 = ri2 with 
In the standard definition of e and 



where uji2 — cji 

ri2 = n(d - ^„)- 

theoretical studies of an oblique impact [14|, n is taken 

at the collision instant, that is, its reorientation during 

the impact is ignored. In experiments, the normal n is 

also determined only once, at the beginning of an impact 
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FIG. 3: Dependence of the angular displacement a = 
arcos [n(0) • n(fc)] of the unit normal n, Eq. 0, on the inci- 
dent angle 7. Open squares and circles are the MD results 
for the model A and B, respectively. Solid and broken lines 
are the corresponding theoretical predictions, a = L" Q{t)dt 
(see text for detail). 



[15J . n = n(0). Neglecting angular velocities (note that 
^1/2(0) = 0) we find for the restitution coefficient: 



g(y-n(0) 
g(0)-n(0) 



Lite) 



Uo) 



d sin ( 



cos a 



Vcos'-f 



(4) 



where we take into account that g(0) • n(0) = — Cti(O) = 
— ycos7 < 0, that ri2(<c) = d and that £,n{tc) < 0. For 
head-on collisions, when 7^-0 and a — > (see Fig. [3]) 
the second term in the r.h.s. of Eq. (j4]) is negligible and 
e is positive. For oblique impacts 7 and a are large and 
the second term prevails, yielding a negative e. 

Hence, the negative restitution coefficient is a conse- 
quence of a significant reorientation of a contact plane 
during a collision. For hard particles with a small col- 
lision duration the reorientation of n is small and may 
be neglected 1J| ; this usually holds true for macroscopic 
bodies. Nanoclusters, however, are very soft particles 
with small Young's modulus [l7[. The duration of their 
impact tc is relatively large and the reorientation of the 
contact plane is significant. 

As it follows from the Eq. ^, the standard definition 
of e characterizes not only the normal motion along n(t) 
(the first term in the r.h.s. of Eq. Q), but also the 
change of the normal n{t) (the second term in the r.h.s. of 
Eq. (HI). Therefore, it is worth to define the restitution 
coefficient, which describes pure normal motion. The 
respective modification of the standard definition reads: 



g(tc) ■ n{tc) 
' g(0) • n(0) 
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(5) 



Here we use Eq. ([3]) for t — tc and take into account that 
h ■ n = for a unit vector n. Note, that the modified 
restitution coefficient e is always positive. Fig. [5J It can 
be also seen from Fig. [5] that the magnitude of e for an 



oblique impact (for large 7) is significantly larger than 
that for a head-on collision. In what follows we explain 
the observed behaviors of e and e using a simple theoret- 
ical model, based on continuum mechanics approach. 

Theory of an oblique impact. Consider a non-inertial 
frame, rotating with the angular velocity H, perpendic- 
ular n, so that n = ft x n. To compute the normal 
force acting between two nanoclusters we apply the im- 
pact theory for macroscopic viscoelastic adhesive spheres 
[j, |5|. It contains the JKR force [16[, which accounts 
for elastic interactions via the Herzian force Fh and for 
adhesive interactions via the Boussinesq force Fb , 



Fh — Fb — 
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It also contains the dissipative force [5|, 



Fd = dri 
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Here, a is the contact radius of the colliding nanoclusters, 
related to the normal displacement ^„ as 



_ 4a2 ISTTcrDa 



(8) 



and D ~ (3/2)(l — v'^jY is the elastic constant with 
the Young modulus Y and the Poisson ratio v. From 
the independent numerical simulations we estimate Y = 
88.3 ela\j and v = 0.396 for inodel A, and Y = 283 GPa 



and V = 0.166 for model B [17|. The surface tension 
a may be expressed via Hamaker constant Ah and the 
equilibrium distance between atoms at the interface zg 
as o- ~ Ah/^Attzq. We obtain a = 0.0246 e/cr|,j and 
0.00289 N/m for the models A and B, respectively. The 
dissipative constant rj, which accounts for the viscoelas- 
ticity of the particles' material [J] is used here as a fitting 
parameter. In the present simulations a good agreement 
is obtained by choosing 77 = 0.95 cr^jy^m/e and 1.62 fs 
for models A and B, respectively. 

In the non-inertial frame, the inertial force must be also 
taken into account. Its normal component reads |18j . 



Fi = 2^vi2 • n{t) - fixi2\n{t)\ 



(9) 



where /i = Nm/2 is the reduced mass of the nanoclusters. 
If wc again neglect the angular velocities of particles in 
the collision (since 1^1 /2{^) = 0), that is, if we assume 
that the two clusters at a contact move together as a 
solid dumbbell, we can exploit the conservation of the 
angular momentum in the form, 



/ir^2^ = ^J'V sin 7 d , 



(10) 



where we take into account that n ■ fi = 0. This yields 
Q,{t) = V sv[\^ d/r'l2{t) and the inertial force. 



hi = — j^ — sm 7 
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Combining Eqs. 
tion for ^„: 



- ((TT|) we obtain the equation of mo- 
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M^Cn + Fh - Fb 



Fd 



{d-U- 



■ sin^ 7 = 0, (12) 



where d' /dt denotes the time derivative in the non- 
inertial frame. Solving Eq. p^ for £,n{t), we obtain e as it 
follows from Eq. ([5]). Taking into account that a = i^{tc) 
we obtain from Eq. Q the relation between the standard 
and modified restitution coefficients, 



e cos a ~ tan 7 sin a 



(13) 



The last equation together with the relation a = 
Jq" n{t)dt may be used to compute the standard coef- 
ficient e. The theoretical predictions for the coefficients 
e and e are shown on the upper and lower panels of Fig. [5] 
respectively. The agreement between our theory, which 
has only one fitting parameter, and MD simulations is 
rather good. We find that the restitution coefficient of 
H-passivated Si nanospheres is well reproduced by our 
macroscopic theory for the incident speed between 20 
m/s and 2405 m/s. If, however, the speed exceeds 2500 
m/s, the nanospheres melt and fuse upon collisions and 
the theory fails to describe the impact. 

We wish to stress that our theoretical model, developed 
for nanoclusters, may be relevant for oblique collisions of 
macroscopic bodies, provided the re-orientation of the 
contact plane during the impact is not negligible. This 
may happen for soft cohesive particles with a low Young 
modulus and large collision time. Relevance of the th eory 
for collisions in wet granular systems is also expected |19| . 

In conclusion, we perform a detailed study of the 
oblique impact of nanoclusters by means of Molecular 
Dynamics and theoretically. In simulations we use two 
models, a simplified one, based on the Lennard- Jones po- 
tential with a cohesive parameter and a realistic model 
for nanoclusters with covalently bonded atoms. We de- 
tect unexpected behavior of the normal restitution coeffi- 
cient e, which becomes negative for large incident angles 
and explain this effect by the reorientation of the contact 
plane in the course of collision. We propose a modified 
definition of the restitution coefficient, e, which describes 
only the normal motion of particles, independently of 
their relative reorientation, and is always positive. A 
simple relation between e and e, that may be helpful 
for experiments is reported. We develop a theoretical 
model for an oblique impact, based on the continuum 
mechanics description of colliding particles, and demon- 
strate that theoretical predictions agree well with simu- 
lation results. Hence, we conclude that the macroscopic 
concepts of elasticity, surface tension and bulk viscosity 
are well applicable for nano-objects of a few hundreds 
atoms. 
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